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ABSTRACT 

We report results from an 18-day ROSAT HRI monitoring campaign on the ultra- 
soft Narrow-Line Seyfert 1 (NLS1) class quasar PHL 1092. This luminous, radio-quiet 
quasar showed strong X-ray variability in a short ROSAT PSPC observation, and 
ROSAT HRI monitoring of the similar object IRAS 13224-3809 revealed extreme 
variability on intermediate timescales. We wanted to determine whether remarkable 
X-ray variability persistently occurs in PHL 1092, and we also wanted to search for 
outstanding variability events that constrain emission processes. Given the large lumi- 
nosity of PHL 1092 (~ 5 x 10 45 erg s _1 in the HRI band), we detect extremely rapid 
and large-amplitude X-ray variability throughout our monitoring. The maximum ob- 
served variability amplitude is a factor of « 14, and in the most rapid variability 
event the HRI count rate increases by a factor of ss 3.8 in a rest-frame time interval of 
< 3580 s. The most rapid event has a rate change of luminosity of > 1.3 x 10 42 erg s~ 2 , 
making it the most extreme such event we are aware of from a radio-quiet quasar. Stan- 
dard 'radiative efficiency limit' arguments imply a radiative efficiency larger than can 
be achieved by accretion onto a Kerr black hole rotating at the maximum plausible 
rate, although we point out that such arguments depend upon the geometry of ini- 
tial radiation release. Relativistic motions of the X-ray source are probably causing 
the radiative efficiency limit to break down; such relativistic motions have also been 
inferred in the similar NLSl-class quasar PKS 0558-504. 
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1 INTRODUCTION 

PHL 1092 (B = 16.7; z = 0.396) is a luminous Narrow-Line 
Seyfert 1 (NLS1) class quasar that is one of the strongest op- 
tical Fe II emitters known (Bergeron & Kunth 1980). Such 
objects have been generally found to have extreme X-ray 
spectral and variability properties, and it appears likely that 
their exceptional X-ray/optical characteristics arise as the 
result of an extreme value of a primary physical parame- 
ter (see Brandt & Boiler 1998 for a recent discussion). This 
primary parameter must ultimately originate from the im- 
mediate vicinity of the supermassive black hole since it can 
strongly influence the energetically-important and rapidly- 
variable X-ray emission. 

The ROSAT PSPC spectrum of PHL 1092 is one of 
the softest ever seen from a quasar (Brandt 1995; Forster 
& Halpern 1996, hereafter FH96; Lawrence et al. 1997). For 
example, a simple power-law fit to the PSPC spectrum gives 



a photon index of F = 4.2 ± 0.5. The poorly-sampled PSPC 
light curve showed remarkably rapid variability for such a 
luminous object. The count rate increased by a factor of ~ 4 
during a 2-day period, and there were weak indications for 
even more rapid variability. No strong spectral variability 
was apparent. FH96 used the 'radiative efficiency limit' of 
Fabian (1979) to argue that the implied radiative efficiency 
was rj > 0.13, and they suggested that a Kerr black hole 
and/or anisotropic emission was implied. The other authors 
cited above interpreted the data in a somewhat more re- 
served manner (see below for discussion), and they found 
r\ ^ 0.02. In this case, a Kerr black hole and/or anisotropic 
emission is not necessarily required. 

We performed an 18-day ROSAT HRI monitoring cam- 
paign on PHL 1092 to further study its X-ray variability 
properties, and here we report the results from our cam- 
paign. Our monitoring goals were (1) to determine whether 
extreme X-ray variability persistently occurs in this quasar 



© 0000 RAS 



2 W.N. Brandt et al. 



m 
-i-i 

a 

o 
u 



o 
u 



C\2 

d 



o 



CO 

I 



1 1 1 1 1 1 — 

PHL 1092 

HRI monitoring 

Linear scale 



n 1 1 r~ 



^ 



u 




mi 



f! 



±jl. 



_i i i i_ 



10 



15 



Logarithmic scale + 4 f 



n 1 1 r - 



♦ ?♦ 



_i i i i_ 



10 

Offset Julian Date 



15 



Figure 1. ROSAT HRI light curve for PHL 1092 obtained during the 18-day monitoring campaign between 1997 July 16 and 1997 Au- 
gust 2. The abscissa label gives the Julian date minus 2450645.12 days. The top panel shows the light curve with a linear ordinate, and 
the bottom panel shows the light curve with a logarithmic ordinate. Each data point is plotted at the middle of the 400 s exposure 
interval from which it was obtained, and the sizes of the exposure intervals lie within the data points themselves. The dashed curve 
indicates the background counting rate within the source extraction circle as a function of time. 



and (2) to search for outstanding variability events even 
more extreme than those seen by the PSPC. As we discuss 
below, goal 2 is important since such variability events can 
place constraints on emission processes and may elucidate 
the origin of the extreme properties of ultrasoft NLS1 more 
generally. This work builds upon our HRI monitoring of the 
similar but lower-luminosity ultrasoft NLS1 IRAS 13224- 
3809 (Boiler et al. 1997, hereafter BBFF). PHL 1092 is ~ 40 
times more luminous than IRAS 13224-3809. 

The Galactic column density towards PHL 1092 is 
(3.6 ±0.2) x 10 20 cm" 2 (Murphy et al. 1996), and the PSPC 
spectrum constrained the intrinsic column density to be 
< 2.0 x 10 20 cm" 2 . We adopt H = 70 km s" 1 Mpc" 1 
and qo = \, and we hence derive a luminosity distance of 
1840 Mpc. When we calculate luminosities below, we shall 
implicitly assume isotropic emission unless stated otherwise. 



2 ROSAT HRI MONITORING RESULTS 

2.1 Observations and spatial analysis 

ROSAT HRI observations of PHL 1092 were performed 
between 1997 July 16 (02:47:15 UT) and 1997 August 2 
(13:53:58 UT). The total exposure time was 109.135 ks, and 
the source was centred on-axis. Up to five separate observa- 
tions were obtained each day, and coverage was reasonably 



good except for a period when PHL 1092 was not observable 
due the position of the Moon (the largest gap in Figure 1). 
All data analysis has been performed using FTOOLS. 

We have added together all the HRI observations to 
produce a master image. The centroid position of PHL 1092 
in the HRI master image is 02000 = 01 h 39 m 56.1 s , #2000 = 
06°19'24.7". This position is in good agreement with the 
precise optical position given by Fanti et al. (1977), and 
there is no significant evidence for X-ray spatial extent after 
errors in attitude correction are taken into consideration (see 
Morse 1994). The source plus background photons at the 
HRI position of PHL 1092 were extracted using a circular 
source cell with a radius of 25 arcsec. The background was 
extracted using an annulus centred on PHL 1092 with an 
inner radius of 40 arcsec and an outer radius of 100 arcsec. 

2.2 Count rate variability 

In Figure 1 we show the full monitoring light curve for 
PHL 1092. Strong X-ray variability is apparent through- 
out the light curve. We have calculated the smallest and 
largest observed HRI count rates to determine the maxi- 
mum variability amplitude. In such calculations, we always 
average over at least 4 of the data points shown in Fig- 
ure 1 to prevent statistical count rate fluctuations from in- 
ducing artificially large variability amplitudes. The smallest 



© 0000 RAS, MNRAS 000, 000-000 



ROSAT monitoring of extreme variability in PHL 1092 3 




2000 



4000 6000 
Time (s) 



8000 



Figure 2. Light curve showing the most rapid observed variabil- 
ity during the HRI monitoring campaign. The HRI count rate for 
PHL 1092 is shown by the large circular data points, and the 
background count rate in the source cell is shown by the small 
square data points. The bin size is 200 s. 



a power-law plus double blackbody model (M2) and a 
power- law plus bremsstrahlung model (M3). Using xspec 
(Arnaud 1996) and pimms (Mukai 1997), we find / M i = 
6.0 x 10 46 erg count" 1 , / M2 = 4.8 x 10 46 erg count" 1 and 
/m3 = 5.2 x 10 46 erg count -1 . In order to be reasonably 
conservative, we shall adopt /m2- 

We find that the 0.2-2.0 keV luminosity of PHL 1092 
varies between (5.4 ± 1.7) x 10 44 erg s" 1 and (7.5 ± 0.8) x 
10 45 erg s _1 during our monitoring observation. If we as- 
sume that the optical continuum is not highly variable on 
timescales of a few weeks (cf. section 3.1 of FH96), the a x 
value for PHL 1092 varies between —1.8 and —1.4 (we cal- 
culate a ox between 3000 A and 2 keV). For comparison, 
the a ox range for optically selected quasars is ~ —1.5 ± 0.1 
(e.g. Laor et al. 1997). For the variability event shown 
in Figure 2, we detect a change in HRI count rate of 
(9.64 ± 2.14) x 10" 2 count s _1 . This corresponds to a lu- 
minosity change of AL = (4.6 ± 1.0) x 10 45 erg s" 1 in a 
rest-frame time interval of At < 3580 s. 



count rate is observed between days 7.2-7.4 in Figure 1, 
and the 6 data points during this span of time give a mean 
count rate of (1.13 ± 0.36) x 10 -2 count s _1 . The largest 
count rate is observed between days 17.0-17.2, and the 4 
data points during this span of time give a mean count rate 
of (1.57±0.16) x 10" 1 count s" 1 . The most probable maxi- 
mum variability amplitude is a factor of 13.9, and maximum 
variability amplitudes in the range 9.5-22.5 are most likely 
(Cauchy distributed; see section 2.3.4 of Eadie et al. 1971). 

In Figure 2 we show the most rapid observed variability 
event. This event occurred around day 8.6 in Figure 1. The 
HRI count rate rose from (3.46 ± 0.98) x 10" 2 count s" 1 
to (1.31 ± 0.19) x 10" 1 count s _1 , thereby increasing by a 
factor of « 3.8 in < 5000 s (< 3580 s in the rest frame of 
PHL 1092). We also detect 6 additional events where the 
HRI count rate exhibits a highly-significant increase or de- 
crease by a factor > 2 in < 1 day. 

2.3 Luminosity variability 

We use the spectrum observed by the ROSAT PSPC to de- 
termine a conversion factor, /, between HRI count rate and 
luminosity. Unlike FH96, we do not simply extrapolate a 
steep (F = 4.2) power-law model down to 0.1 keV to de- 
rive a 0.1-2.0 keV luminosity. The extrapolation of a steep 
power-law spectrum to low energies where Galactic absorp- 
tion is important can lead to unphysically large, or at least 
highly uncertain, luminosities. For such models most of the 
luminosity lies between 0.1-0.2 keV where few counts are 
actually detected.^] 

We have instead considered the luminosity from 0.2- 
2.0 keV in the observed frame since, for the Galactic col- 
umn density relevant here, this quantity is much less de- 
pendent upon the details of the spectral model. We have 
considered three spectral models that give statistically ac- 
ceptable fits to the PSPC data: a power-law model (Ml), 

* This is the cause of the much higher radiative efficiency de- 
rived by FH96 as compared to Brandt (1995) and Lawrence et al. 
(1997). 



3 DISCUSSION 

3.1 Radiative efficiency limit arguments 

The variability event shown in Figure 2 has ^ > 1.3 x 
10 42 erg s -2 , making it the most extreme such event we are 
aware of from a radio-quiet quasar (compare with Remil- 
lard et al. 1991 and table 3 of FH96). This can be further 
quantified by employing the radiative efficiency (77) limit: 
77 > 4.8 x 10~ 43 ^H (Fabian 1979). Straightforward applica- 
tion of the limit gives an extremely high r\ > 0.62 ± 0.13. 
For comparison, optimal accretion onto a Kerr black hole 
rotating at the maximum plausible rate gives only r/ ~ 0.3 
(see Thorne 1974). 

The extreme radiative efficiency derived above provides 
motivation to critically examine the assumptions underly- 
ing the radiative efficiency limit. The standard derivation 
assumes that the radiation release associated with an ob- 
served luminosity outburst occurs entirely at the centre of 
the emission region. If the radiation release is more uniform, 
the rapid emission of photons from the outer few Thom- 
son depths facing the observer can invalidate the standard 
derivation (see Appendix A). A rapid rise in the flux from a 
source, such as we have for PHL 1092 (and such as has been 
typically used in the literature), need not therefore involve 
the whole source and no limit applies. 

We can, however, recover the situation if either (a) we 
assume restrictions upon the manner of radiation release or 
(b) we can place an upper limit upon the Thomson depth of 
the emission region. For case (a), we return to the original 
efficiency limit if all parts of the region emitting the sharp 
flux increase are in causal contact with each other (rather 
than with some remote centre of an idealized sphere; see Ap- 
pendix A). For case (b), we have from Appendix A when tt 
is large that At = #^ = $7^; where $ is a geometrical fac- 
tor of order unity. Thus tt = tf^j • H we can place an upper 
limit on R we also place an upper limit on tt. The emission 
dominating the ROSAT band is from the soft X-ray excess, 
and this emission is thought to be associated with the inner 
accretion disk. If we first make the assumption of blackbody 
emission we obtain R < 4 x W 11 L^T' 2 cm where L44 is 
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the total blackbody luminosity in units of 10 erg s _1 and 
T% is the blackbody temperature in units of 10 6 K. We have 
1/44 ~ 50, and Te is unlikely to be smaller than ~ 0.35 based 
on accretion disk theory and observations of the soft X-ray 
excess in similar objects (see equation 3.20 of Peterson 1997; 
measured soft X-ray excess temperatures are usually larger 
than this but will be strongly affected by relativistic effects 
if the black hole is Kerr). We then find R < 2.3 x 10 13 cm 
and tt ~ 0.5$. While the derived R constraint is physically 
plausible (of order a couple gravitational radii for a black 
hole mass of ps 10 8 Mq), the low Thomson depth is incon- 
sistent with the assumption that tt is large. Thus when the 
source size is determined by blackbody radiation, one can- 
not escape the efficiency limit by invoking large tt as is done 
in Appendix A. Additional opacity sources, such as free-free 
absorption, are expected in the case of blackbody or quasi- 
blackbody emission. These have the effect of strengthening 
the above argument (tt is replaced by an effective Trot > tt 
that includes the additional opacity sources). 

Another possibility for case (b) is that the emission is 
due to Comptonization (which is plausibly the fastest emis- 
sion process in this band). In this case, a cloud of gas bathed 
in soft photons is impulsively heated. The soft photons must 
originate at energies below the HRI band for large variabil- 
ity to be seen. Large spectral variations will occur as the 
outgoing photons are successively Compton upscattered be- 
fore escaping unless the electron temperature is large and 
tt is small (see Guilbert, Fabian & Ross 1982). The relative 
change in photon energy per scattering is 4kT/mc 2 , and 
therefore kT ~ mc 2 /4 in order that many of the first burst 
of photons, which scatter at most once, appear in the HRI 
band. Such a high temperature requires that the Thomson 
depth of the whole source is small. The equilibrium Compton 
scattering formulae given by Zdziarski (1985; e.g. his equa- 
tion 5; we correct his equation 6) require that tt ~ 0.04 if 
kT/mc 2 ~ 0.25 in order to obtain the steep mean photon 
index from the source of T = 4.2 ± 0.5. Again this upper 
limit on tt allows a lower bound to be placed upon 77. 

The precise value of the coefficient in the efficiency limit 
does depend upon geometry and can be uncertain by a fac- 
tor of a few. We have seen that plausible constraints on 
radiation mechanisms make it unlikely that the much larger 
uncertainty theoretically possible for a high tt sphere is re- 
alised. We note further that there are likely to be inefficien- 
cies in the conversion of the accretion energy into radiation, 
with some energy ending in kinetic or magnetic form. We 
conclude that the rapid variation seen in Figure 2 requires 
an efficiency of at least « 0.10 and more probably ~ 0.60. 

3.2 Relativistic X-ray boosting in ultrasoft NLS1? 

A straightforward application of the radiative efficiency limit 
to our most rapid observed variability implies r\ > 0.62 ± 
0.13. This appears to be substantially larger than can be 
explained without allowing for relativistic effects, but mod- 
erate relativistic boosting can easily explain our data (see 
Guilbert, Fabian & Rees 1983 and BBFF). The relativistic 
bulk motions required are of order 0.3c and might plausibly 
arise in the inner accretion disk. Relativistic X-ray boost- 
ing associated with a strong jet appears less likely given the 
spectral character of the soft X-ray emission and the fact 
that PHL 1092 is radio quiet. 



Evidence for relativistic X-ray variability enhance- 
ment has been found in two other NLSl-class objects to 
date: PKS 0558-504 and IRAS 13224-3809. PKS 0558-504 
showed a rapid X-ray flare that implied relativistic motions 
(Remillard et al. 1991), while IRAS 13224-3809 shows per- 
sistent, giant-amplitude, nonlinear variability that is most 
naturally explained via relativistic effects (BBFF). Our data 
on PHL 1092 now add further weight to the idea that un- 
usually strong relativistic effects may be present in many 
ultrasoft NLS1. This hypothesis can be further examined 
by searching for X-ray spectral changes during putative rel- 
ativistic variability events. XMM has the large collecting 
area and low-energy sensitivity needed for this work, and 
we hope to perform such observations. 
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Figure Al. The photon flux detected at an observer versus 
time t for 10 logarithmically spaced values of tt from Tp = 0.1 to 
"Pr = 40. Smaller Thomson depths correspond to the curves with 
lower maximum values. 



APPENDIX A: DEPENDENCE OF THE 
RADIATIVE EFFICIENCY LIMIT UPON THE 
GEOMETRY OF RADIATION RELEASE 

The radiative efficiency limit was derived by arguing that a 
radiatively-inefiicient source must have such a high particle 
density that a rapid release of radiation in the source appears 
to an observer to be slowed down by Thomson scattering. 
For a luminosity outburst of amplitude AL with a rise time 
of At, the radiative efficiency r\ is defined by the equation 



ALAt = rjMc 



(Al) 



where M is the mass involved in the outburst and c is the 
speed of light. M w m p nV where m p is the mass of the 
proton, n is the proton number density (also assumed to be 
the electron number density), and V is the volume of the 
emission region. The standard derivation of the radiative 
efficiency limit (Fabian 1979) assumes a uniform, spherical 
emission region with the release of radiation localized at its 
centre, and it also assumes that relativistic Doppler boosting 
and light bending are unimportant. If the emission region 
has a significant Thomson depth (tt > 1), At must satisfy 



At > (1 + 7T)-, 
c 



(A2) 



where <tt is the Thomson cross section and R is the radius 
of the region. The limit on r\ arises due to the competition 
between the light crossing time and the photon diffusion 
time. By combining equations (Al) and (A2) one obtains 



>» AL tl \ 
^ ~~At 



/(tt) oc 



(1 + tt) 

TT 



(A3) 



/(tt) has the asymptotic behaviour 



/(tt) oc 



TT 



for tt <C 1 (light crossing time dominates) 
for tt 3> 1 (photon diffusion time dominates). 



This means that there exists a minimum value of /(tt) (4 
when tt = 1) and hence of r\ for a given AL and At. 

In order to relax some of the assumptions entering the 
standard analytic derivation, we simulated photon diffusion 
from 'clouds' of moderate and high Thomson depths using 



Monte Carlo techniques. We first considered photon diffu- 
sion for the case described above with high Thomson depth 
and found our results to be consistent with the analytic re- 
sults published by Sunyaev & Titarchuk (1980). We then 
considered a simple model of a uniform, spherical cloud but 
allowed radiation to be produced instantaneously through- 
out the cloud in a uniform manner. The results of the sim- 
ulations for different Thomson depths are presented in Fig- 
ure Al. The curves show the time dependence of the ob- 
served photon flux ^ oc L for a wide range of tt- There 
is an important difference in the behaviour of these curves 
as compared to the case when the radiation release occurs 
entirely at the centre of the cloud: when tt becomes large 
there is now no reduction in the rate of increase of -rr at 

at 

the start of the outburst. In other words, at the start of 
the outburst oc ^ increases monotonically as tt is in- 
creased (even for t t 1). This result can be understood by 
realizing that the radiation observed near the start of the 
outburst comes mostly from the outer few Thomson depths 
of the cloud facing the observer. The relevant 'cup-shaped' 
emission region at the start of the outburst is constrained 
by the surface of the cloud and the surface of equal photon 
arrival time. An upper limit to its volume for the diffusion 
dominated case is Ku P < irR 2 h, where h w Appropriate 
modification of equation (A3) leads to 



where 

/(t t ) oc t t 2 



(A4) 



for tt ^> 1 (photon diffusion time dominates). 



The lower limit on the radiative efficiency in this case is a 
trivial r\ > 0. 

The assumption that radiation is released instanta- 
neously throughout the cloud is obviously an unrealistic 
one because it violates causality. To investigate the effect of 
causality, we considered a 'trigger signal' propagating out- 
wards from the centre of the cloud with speed c. Radiation 
is released in a volume element of the cloud only after the 
trigger signal reaches it. Inclusion of this effect in the sim- 
ulations does not change the general behaviour described 
above; at the start of the outburst, 



777^ 



increases monoton- 



ically as tt is increased. Again the only constraint on the 
radiative efficiency is the trivial r\ > 0. 

A nonzero lower limit for r\ can be derived if one consid- 
ers somewhat different definitions for AL and At. For exam- 
ple, one could define the characteristic variation timescale 
At as the time (measured from the observed start of the 
outburst) it takes for ~ 90 per cent of the outburst photons 
reach the observer. Similarly AL could be defined as 90 per 
cent of the energy liberated in the outburst divided by At. 
There is a nonzero lower limit for r/ with the above defini- 
tions (this has been checked numerically). However, in order 
to reliably obtain an efficiency bound in this way one must 
be able to determine the overall profile of an outburst. In 
particular, one must be able to measure any decaying 'tail' 
with reasonably high precision. This is usually not possible 
due to observational constraints and the complexity of ac- 
tive galaxy light curves (e.g. outbursts often overlap each 
other). 

In summary, the radiative efficiency limit is quite sen- 
sitive to the geometry of radiation release within the emis- 
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sion region. Without constraints on this geometry or the 
Thomson depth of the emission region, it is difficult to place 
meaningful constraints on r\. 
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